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ABSTRACT 

We present the results of a spectroscopic monitoring campaign devoted to three O-type stars 
in the Carina nebula. We derive the full SB2 orbital solution of the binary system Tr 16- 
1 12, an exceptional dissymmetrical system consisting of an 05.5-6 V((f+?p)) primary and a 
B2 V-III secondary. We also report on low-amplitude brightness variations in Tr 16-112 that 
are likely due to the ellipsoidal shape of the 05.5-6 primary revolving in an eccentric orbit 
around the system's centre of mass. We detect for the first time a clear SB2 binary signature 
in the spectrum of HD 93343 (08 + 08), although our data are not sufficient to establish an 
orbital solution. This system also displays low amplitude photometric modulations. On the 
other hand, no indication of multiplicity is found in the optical spectra of HD 93250. Finally, 
we discuss the general properties of multiple massive stars in the Carina OBI association. 

Key words: stars: individual: Tr 16-112 - stars: individual: HD 93343 - stars: individual: 
HD 93250 - stars: binaries: spectroscopic - stars: early-type - stars: fundamental parameters 



1 INTRODUCTION 

The investigation of the multiplicity of early-type stars in very 
young open clusters is a powerful tool to better understand many 
problems related to massive star research. First of all, such studies 
allow us to enlarge the - still rather small - sample of massive stars 
with well determined fundamental properties such as masses and 
radii. Moreover, the statistical properties of the binary systems also 
constrain the formation scenarios of massive stars (see e.g. |Sana"etl 
|al.|2008| i. 

The region around 77 Carinae harbours several very young 
open clusters and is therefore extremely rich in early-type stars. 
As a consequence, this area has been the subject of many investi- 
gations over the last thirty years. In this context and in the frame- 
work of the so-called X-Mega project ( [Corcoran et al.|1999[ l, aim- 
ing (among other things) at an intensive monitoring of the region 
around 77 Carinae with the ROSAT X-ray satellite, a large multi- 
wavelength observing effort was devoted to the Carina Nebula and 



* Based on observations collected at the European Southern Observatory 
(La Silla, Chile), at Complejo Astronomico El Leoncito (Argentina), at the 
Cerro Tololo Inter-American Observatory (CTIO) and with XMM-Newton, 
an ESA Science Mission with instruments and contributions directly funded 
by ESA Member States and the USA (NASA), 
t Research Associate FRS/FNRS (Belgium) 
X Postdoctoral Researcher FRS/FNRS (Belgium) 
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the Trumpler 16 cluster in particular. One of the prime objectives 
was to obtain accurate ephemerides of coUiding-wind early-type 
binary systems to interpret their X-ray light curves. The spectro- 
scopic monitoring of these objects started in the late 1990's and led 



to a series of p apers ( ^ Albacete Colombo et al. 200I| |2002, ,Mor- 



Irell et al.|200Tl|Rauw et al.|2001||Naze et al.|2005]|Niemela et al. 
|2006[ l reporting the results of dedicated investigations of the prop- 
erties of a number of O-type binaries in the Carina Nebula. This 
paper is the seventh in this series and is devoted to three objects: 
the binary systems Tr 16-1 12 and HD 93343 as well as the likely 
single star HD 93250. Furthermore, we briefly discuss the general 
properties of O-type binaries in and around Trumpler 16 as derived 
in the framework of the present campaign and other related studies. 



2 OBSERVATIONS 

2.1 Medium resolution spectroscopy 

Six medium resolution spectra of Tr 16-112 and of HD 93343 in 
the wavelength range 3850 — 4800 A were gathered in 1997 with 
the ESO 1.5 m telescope equipped with a Boiler & Chivens (B&C) 
Cassegrain spectrograph (see Rauw et al.|2001| for a detailed de- 
scription of the instrumentation and the data reduction). 

Several spectra of Tr 16-1 12 were obtained in 1986 with the 
Cassegrain spectrograph at the CTIO 1 m telescope equipped with a 
two-dimensional, photon-counting detector (2D-frutti). The wave- 
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length coverage was from 3800 to 5000 A with a 3 pixel resolu- 
tion of 1.5 A. The S/N ratios were about 50. Several low-resolution 
spectra of the same star were obtained in 1993, 2007 and 2008 
at Complejo Astronomico El Leoncito (CasleoQ Argentina. The 
spectra were taken with the modified REOSC SElj^ spectrograph 
in simple dispersion Cassegrain mode attached to the 2. 15 m Jorge 
Sahade telescope. The dispersion was 1 .64 A per pixel. The data 
cover the wavelength range from 3800 to 4800 A and have a S/N of 
about 200. The spectra were reduced at Facultad de Ciencias As- 
tronomicas y Geoffsicas (Universidad Nacional de La Plata) with 
the IRAljf] software version 2.10. Radial velocities of the Hell 
A 4542 and A 4686 lines were obtained by fitting a Gaussian profile 
to the observed lines. 



2.2 High resolution echelle spectroscopy 

During several observing campaigns (May 1999, May 2000, May 
2001, March & April 2002 and May 2004), we took a series of 
echelle spectra of our targets (27 for Tr 16-112, 9 for HD 93343 
and 6 for HD 93250) with the Fiber-fed Extended Range Optical 
Spectrograph (FEROS) mounted first on the ESO 1.5 m (until Oc- 
tober 2002) and later on the ESO/MPG 2.2 m telescopes at La Silla. 
These spectra cover the wavelength domain from about 3750 to 
9000 A with a resolving power of 48000. The detector was an EEV 
CCD with 2048 x 4096 pixels of size 15 x 15 fim'^. The exposure 
times were typically 30 min and the typical S/N ratios were about 
150. The data were reduced using the specific context under the MI- 
DAS environment along with an improved reduction pipeline (see 
[Sana et al.|2003^ . Finally, the FEROS spectra were normalized to 
the continuum using polynomials. 

High-resolution echelle spectra were also obtained at Casleo 
in 1994 and 2008 with the REOSC SEL instrument in echelle 
mode attached to the 2.15 m telescope. The detector was a Tek 
1024 X 1024 pixel CCD. The spectra cover the wavelength range 
from 3600 to 6000 A with a spectral resolving power of 18000. 
Typical S/N ratios are about 80. The REOSC echelle spectra 
were reduced using IRAF routines. Individual echelle orders were 
normalized using carefully chosen continuum windows. 



2.3 Photometry of Tr 16-112 

CCD photometry of Tr 16-1 12 was performed during the nights of 
February 21-27, 2007, using the CH-250 camera with a PM512 
CCD chip, attached to the 0.6 m Helen Sawyer Hogg telescope 
(f/15 Cassegrain) at Casleo. We used a standard Johnson-Cousins 
B VR I filter set for the observations, although the bulk of the data 
were taken in the V band. The CCD images were calibrated in the 
standard way using bias, dark and flat-field frames and finally the 
photometry was derived using the APPHOT/IRAF package. The 
typical uncertainty on the differential l/-band photometry of Tr 16- 
112 is about 0.005 mag. 



^ Casleo is operated under agreement between CONICET and the National 
Universities of La Plata, Cordoba and San Juan. 

^ Spectrograph Echelle Liege (jointly built by REOSC and Liege Observa- 
tory and on long term loan from the latter). 

^ IRAF is distributed by NOAO, operated by AURA, Inc., under agreement 
with NSR 



3 INTERSTELLAR LINES 

We added up all the FEROS spectra of our targets in the wave- 
length domain 3870 - 3950 A to derive high quality spectra of the 
interstellar He I A 3889 and Call K lines. For all three stars, the 
profiles of the Call K lines look remarkably similar to the plots 
given by |Walbornlj 19821 . The high signal to noise ratio of our mean 
spectra allows us to further detect some weak (equivalent widths of 
a few mA) higher velocity components that were not seen in the 
older data. In the case of Tr 16-112, we definitely find a compo- 
nent at — 143kms~^ and possibly another one at +205 kms^^. In 
HD 93343, we clearly identify an absorption at — 75kms~^, whilst 
in HD 93250 we detect a component at — 98kms~^. All three stars 
display a moderately prominent interstellar He I A 3889 absorption 
line at a heliocentric velocity of — 28.0kms~^ (for Tr 16-112 and 
HD 93250) and -25.4 km s"^ (HD 93343). The equivalent widths 
of these features are 0.010 (HD 93343), 0.045 (HD 93250) and 
0.080A(Tr 16-112). 



4 THE O + B BINARY TR 16-112 

Trl6-112 {= CPD-59°2641 = CD -59° 3310, |Feinstein, MaF^ 
|raco & Muzzio|1973^ is a member of the Tr 16 cluster. Based on 
10 low-resolution spectra, fLevato et al.| ( |199T| l classified Trl6- 
112 as an SBl binary system with an orbital period of 4.02 days. 
Their preliminary orbital solution yielded a significant eccentricity 
of 0.23. |Luna et al.|(2003^ subsequently obtained 5 additional RV 
data points from echelle spectra and re-computed an SBl orbital 
solution with a period of 4.080 days and e — 0.25 ± 0.03. 

4.1 The spectral types 

We have applied a disentangling method based on the principles 
discussed by |Gonzalez & LevatoH20()6l l to the normalized FEROS 
spectra of Tr 16-1 12. For this purpose, we have selected two wave- 
length domains within the spectral range covered by the FEROS 
data: our violet and blue domains cover respectively the wavelength 
intervals [4000, 4250] and [4410,4790]. Disentangling the spectra 
of the components of this star was a challenging task that pushed 
our code to its limits. Indeed, the secondary star contributes only 
about 8 percent of the total light of the system making it rather 
difficult to accurately follow its spectral signature. 

The resulting separated spectra are displayed in Fig.[T| Al- 
though the disentangled secondary spectrum is admittedly of poor 
quality, these spectra reveal some interesting features. The appar- 
ent equivalent widths (i.e. with respect to the combined primary -l- 
secondary continuum) of the He I A 447 1 line as determined from 
the mean disentangled spectra are 0.325 ± 0.005 and 0.073 A for 
the primary and secondary respectively. For the He II A 4542 and 
A 4686 lines, the apparent primary's equivalent widths (EWs) are 
0.668 ± 0.005 and 0.426 A respectively. None of the latter lines is 
seen in the secondary spectrum (see Fig.[T](. 

The He I A 447 1/He II A 4542 EW ratio (logW = -0.313 ± 
0.007, see |Conti|1973b| and [Mathys|1988| l places the primary star at 
the border between spectral types 05.5 and 06. The disentangled 
primary spectrum reveals a weak, but definite, emission blend of 
NlII AA 4634-41 as well as a comparably strong Cm AA 4647-50 
emission and a weak CiV A 4658 P-Cygni profil^ The fact that 
the strength of the Cm AA 4647-50 and NlII AA 4634-41 blends 

^ From its wavelength, we can actually rule out that the absorption compo- 
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Figure 1. The normalized disentangled violet (left) and blue (right) spectra of the primary and secondary components of Tr 16-1 12. The disentangled spectra 
were normalized accounting for a brightness ratio (secondary/primary) of 0.09. For clarity, the primary spectrum was shifted upwards by 0.2 continuum units. 



are very much comparable would put the primary into the scarce 
category of Of?p stars (see |Walborri[T972l[l^73l and |Naze, Walboml 
|& Martins, 2008 ) . We caution however that the known members of 
the Of?p category exhibit supergiant-like features and often dis- 
play strong spectroscopic variability (see |Naze et al.|2008 |. In our 
data, we find no evidence for such features in the spectrum of the 
primary star of Tr 16-112. The Hell A4686 absorption displays a 
slightly asymmetric profile with some hints of a weak P-Cygni type 
emission component. Some evidence for weak S IV AA 4486, 4504 
emissions is found as well as a very weak Si IV A 41 16 emission, 
whilst Si IV A 4089 is clearly seen in absorption. These features mo- 
tivate a spectral classification of the primary as 05.5-6 V((f'"?p)). 

Using the Fourier transform method (see |Sim6n-Diaz & Her-| 
|rero||2007] l, we have determined the projectetl rotational velocity 
V smi = (170 ± 5) km s"^ from the He II AA 4200, 4542 and O III 
A 5592 lines in the disentangled primary spectrum. This value is 
slightly lower than the values previously determined from th e com- 
posite lUE spectra of the system (195 and 184kms~^, see Penny 



Table 1. Spectral lines used for the radial velocity determination within the 
disentangling procedure. 



|1996| and |Howarth et al.|fl997| respectively). One reason for this 
difference could be blending effects with the secondary's spectrum 
in the lUE data, as Penny mentioned the possible presence of the 
secondary signature in the cross-correlation function. 

The lack of He II lines in the spectrum of the secondary in- 
dicates that this is a B-type star. The most prominent features in 
the disentangled secondary spectrum are the He I lines (AA 4009, 
4026, 4121, 4144, 4471, 4713) as well as US. Somewhat weaker 
features are due to the C III -H O II blends at 4070 and 4650 A, as 
well as the Si III AA 4552, 4568, Mg II A 448 1 , N II A 4530 and O II 
A 41 10 lines. The presence of a very weak Si IV A 4089 absorption 
cannot be ruled out. The relative strength of the lines of the sili- 
con ions and the strength of the Mg II line compared to He I A 447 1 
suggest a spectral type B2 (or slightly earlier) with an uncertainty 
of about half a subclass (see |Walbom & Fitzpatrick|1990[ l. The rel- 
ative strength of the N II A 4530 line argues for a V-III luminosity 
class (see the luminosity effects at spectral type B2 in |Walbom &| 
|Fitzpatrick|i990^ . 

To evaluate the brightness ratio between the components of 



nent of this P-Cygni profile would be due to a blend with the Si IV A 4654 
absorption. 
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the system, we have compared the equivalent widths of prominent 
lines as measured in the combined spectrum to the typical EWs of 
presumably single stars of the same spectral type. For the primary, 
we used the Hel AA4026, 4471 as well as Hell A A 4200, 4542 
lines and we took the typical values of the EWs from the paper of 
|Contj]|T973a|. For the secondary, our comparison was based on the 
He I AA 4026, 447 1 and Mg II A 448 1 lines and the typical values of 
the EWs of B2 V stars given by Didelon ( 1982). This comparison 
shows that the secondary contributes about 8 ± 3 percent of the 
total light of Tr 16-1 12 in the blue- violet spectral domain (i.e. the 
brightness ratio secondary/primary amounts to 0.09). 

The disentangling code offers the possibility of simultane- 
ously determining the radial velocities (RVs) of the binary com- 
ponents. For this purpose, we used the lines listed in Table[T] 

4.2 The orbital solution 

Since the spectral signature of the secondary component of Tr 16- 
112 is rather weak, we started by treating the system as an SBl 
binary. We thus measured the radial velocities (RVs) of the pri- 
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Figure 2. Top panel: periodogram (semi-amplitude spectrum) of the entire 
time series of 98 radial velocities of the primary component of Tr 16-112 
computed with the method of |Heck et al.|{T985^ . Bottom panel: spectral 
window corresponding to our time series. 



mary component by fitting gaussian line profiles to tlie strongest 
lines. It was only on a few FEROS spectra tiiat we directly noted 
the presence of tiie secondary's spectral signature in the He I hnes. 
In the latter cases, we simultaneously fitted two gaussian profiles. 
Since the secondary is of a significantly later spectral type than 
the primary, the He II lines were not affected and we thus consider 
that the radial velocities of the He II lines best reflect the motion of 
the primary of Tr 16-112. In this way, we obtained a time series of 
48 RV measurements from high-resolution echelle spectra, span- 
ning almost 9 years (3259 days). Note that the RV dispersion of the 
diffuse intersteUar band (DIB) at 5780 A in the FEROS spectra of 
Tr 16-112, was found to be 2.3 kms'^. This value can be seen as 
an estimate of the uncertainty on the RV of a constant spectral line 
in the spectrum of this star. 

We then used the generalized periodogram Fourier technique 



discussed by ^Heck, Manfroid & Mersch^ (1985) and Gosset et al.| 
( |2001^ as well as the trial period method of |Lafler & Kinman ^1965^ 



to perform a period search on the primary radial velocities. Previ- 
ous attempts to determine the orbital period from a subset of the 
data discussed here revealed a severe ambiguity between a period 
of almost exactly 4.0 days and its alias at 1.3 days. To get rid of 
the aliasing problem, we thus obtained several spectra during the 
same night in the course of the May 2001 FEROS campaign and 
during the Casleo campaigns in 2008. The periodogram built from 
the new FEROS -I- Casleo data now clearly reveals its highest peak 
at a frequency of (0.24902 ± 0.00003) d"^ corresponding to an 
orbital period of (4.0157 ± 0.0005) days. Our value is in excellent 
agreement with the 4.02 day period originally proposed by [Levatol 
letal. (19911. 



Table 2. Radial velocities of Tr 16-1 12 as measured on our high-resolution 
echelle spectra. The first column yields the date of the observation. The sec- 
ond column indicates the primary RV as inferred from the direct measure- 
ment of the He II lines. Columns 3 and 4 yield the primary and secondary 
RVs as inferred from the disentangling method. Finally, the last column in- 
dicates the instrument that was used for the observation: FEROS (F) or the 
Casleo (C) SEL echelle spectrograph. 
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To check this result, we added a set of 56 older, lower resolu- 
tion, RV measurements to our high-resolution data, thereby extend- 
ing the time span of our data set to 24 years (8783 days). The anal- 
yses of this extended time series clearly confirm the conclusions 



from our high-resolution data set (see Fig.|2]), slightly reducing the 
uncertainty on the period to (4.0157 ± 0.0002) days. 

As a next step, we used the Liege Orbital Solution Package 
[Wolfe, Horak & Storer|p967) re- 



(LOSP) based on the method of 
vised by [Sana, Gosset & Rauw j2006[ >, to compute the SB 1 and 
SB2 orbital solutions of Tr 16-1 12. For the SBl solution, we used 
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Table 3. Orbital parameters of the SBl (left) and SB2 (right) solutions of Tr 16-1 12. To yields the time of periastron passage. The quoted uncertainties 
correspond to the 1 a error bars. 
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Figure 3. Left: SBl radial velocity curve of Tr 16-112 as derived from the RVs of the Hell lines in the primary spectrum. The filled circles stand for the 
FEROS data, whilst the open triangles indicate RVs derived from the Casleo echelle spectra. Right: SB2 radial velocity curve of Tr 16-1 12 as derived from the 
RVs obtained through the disentangling of the FEROS spectra. Filled symbols stand for the primary whilst the open symbols indicate the secondary star. The 
orbital solutions shown are described in Table[3] Phase (j> = 0.0 corresponds to the periastron passage. 



the RVs of the primary's He II lines as measured on the entire set 
of high-resolution echelle spectra, whilst we adopted the RVs de- 
rived from the disentangling method applied to the FEROS spectra 
for the SB2 solution. The corresponding RV curves are shown in 
Fig.lSland the parameters of the orbital solutions are listed in Ta- 
blelsTNote that leaving the orbital period as a free parameter did 
not change the quality of the orbital solution. 



The mass ratio (primary/secondary) is found to be 3.49±0.08. 
To our knowledge, this is one of the largest mass ratios found in 
an OB binary system to date. Our orbital solution confirms the 
existence of a significant eccentricity, although it is less extreme 
than previously thought. Note that most parameters agree reason- 
ably well between the SBl and SB2 orbital solutions, except for 
the longitude of periastron which has a rather large uncertainty. 
Comparing the relatively large minimum mass of the secondary 
with the typical masses of B2 stars suggests that the orbital inclina- 
tion could be as large as 60°. Under these circumstances, grazing 
eclipses could a priori be possible. 



4.3 Photometric data 

Our photometric data reveal that Tr 16-1 12 displays low level vari- 
ations in its light curve with an amplitude of about 0.03 mag. When 
folded with the ephemeris found from the radial velocity curve, 
we note that our observations actually miss the phase when sec- 
ondary eclipse would be expected (around (j> = 0.842) but cover 
the phase where primary eclipse would occur (0 = 0.251). The 
photometric data reveal no trace of an eclipse at this orbital phase 
and we thus conclude that Tr 16-112 does not display photometric 
eclipses. Rather the low-level photometric variability must be due 
to the ellipsoidal shapes of the stars. 

We analysed the light curve of this system by means of the 
NIGHTFALL cod^developed and maintained by R. Wichmann, M. 
Kuster and R Risse. This software uses a Wilson-Devinney like it- 
erative approach to solve the light curves of close binary systems in 



http://www.hs.uni-hamburg.de/DEAns/perAVichmann/ 
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Table 4. Parameters of the best fit model for the light curve of Tr 16-112. 
The error bars correspond to the 90% confidence intervals. The average 
residual of the fit amounts to |0 — C| = 0.005 mag. 
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the Roche geometry, generahzed to deal with the more complex sit- 
uation of eccentric binaries where the sizes of the Roche lobes (and 
hence the values of the filling factors) change with orbital phase 
(see |Kopal|1989i ). Despite the limited phase coverage and the small 
amplitude of the variations, we can attempt to get a decent fit of 
the data by varying the relevant model parameters which are in this 
case the orbital inclination, as well as the primary and secondary 
filling factors, defined as the ratio between the polar radius of the 
star and the polar radius of the Roche lobe at periastron passage. 

We set the effective temperature of the primary star to 39000 K 
( [Martins, Schaerer & Hillier|2005a^ whilst the secondary temper- 
ature was fixed to 22000 K according to the |Schrnid t-Kaler ( 1982^ 
calibration, which seems reasonable for B2 main-sequence stars 
P<itzpatrick & Massa|2005j . The mass ratio, the orbital eccentric- 
ity as well as the longitude of periastron were fixed at the values 
derived from the radial velocity curve (see Table[3]l. The NIGHT- 
FALL code was applied to the data in the V-band since they are 
more numerous (996 data points) and of the highest quality. Due 
to the limited phase coverage of the data and the small range of 
variability, the parameters are actually not that well constrained. To 
evaluate the uncertainties on the model parameters, we systemati- 
cally explored the parameter space, fixing one parameter at once, 
and searching for the combinations of the other two parameters that 
yield the lowest residuals. 

The best fit value of the orbital inclination is found to be 54°. 
This then corresponds to absolute masses of (41.0 ± 2.5) M© and 
(11.7 ± 0.6) M© for the primary and secondary respectively. 

Given the rather large difference in temperature (and thus in 
surface brightness) between the primary and the secondary, the fit 
of the light curve is not very sensitive to the value of the secondary 
filling factor. As a result, the light curve does not allow us to effec- 
tively constrain the brightness ratio of the two stars independently 
of the results obtained from the spectroscopy. The situation is better 
for the primary (which dominates the light of the system). For this 
component, we can constrain the filling factor to be in the range 
0.72 to 0.74 with a preferred value of 0.73. 

We emphasize that the parameters in Table|4] should be taken 
as a first estimate only. A more extensive, higher quality (ideally 
space-borne), photometric campaign would be needed to fill in the 
gaps in the light curve and ascertain these numbers. Here, we only 
wish to stress that the observed light curve can be explained through 
ellipsoidal variations. Pending confirmation of these parameters by 
future studies, we note that the primary is slightly more massive, 
significantly larger and brighter than a 'typical' main sequence star 
of the same spectral type. In fact, these parameters would be more 
representative of luminosity class III objects ^Martins et al.|2005al l, 
although this would be at odds with the presence of the strong He II 
A 4686 absorption in the primary spectrum. This would then sug- 
gest that Tr 16-1 12 is actually a slightly evolved object, contrary to 




0.5 1 

0(P = 4.01 57 days) 

Figure 4. Light curve of Tr 16-1 12 in the V-band (top panel). The solid 
curve yields the best fit model where the variations are due to the ellipsoidal 
shape of the primary star filling part of its critical volume in an eccentric bi- 
nary. The parameters of the fit are listed in Table[4] The residuals are shown 
in the bottom panel. 



systems such as Tr 16-104 ( |Rauw et al.|200U or FO 15 jNiemela"etl 
al. 2006 1. Another point worth mentioning is that if the above pa- 
rameters are indeed confirmed and if the reddening law towards this 
star is characterized by the typical value of the Trumpler 16 cluster 
(Rv = 3.48 ± 0.33, see |Carraro et al.||2004| l, then the distance 
modulus of Tr 16-1 12 would be 13.0 ± 0.3 (i.e. corresponding to 
a distance of 4.0 kpc), almost twice the distance inferred from the 
study of several eclipsing binaries in Tr 16 (e.g. Tr 16-l, |Freyham-| 
|mer et al.|20"0T] and Tr 16-104, |Rauw et al.|200l |l. Whilst this is a 
surprising result, we have to recall that our line of sight towards the 
Carina region is actually roughly along the Sagittarius-Carina spiral 
arm, so that several distinct populations of massive stars could be 
seen projected on the same area of the sky. In the case of the eclips- 
ing binary FO 15 jNiemela et al.|2006] l, a similar contradiction was 
solved by showing that the latter star is subject to a highly anoma- 
lous extinction with Rv = 4.15. We note that for Tr 16-1 12, there 
is currently no evidence for such a large value of Rv ■ In addition, 
since the E{B ~V) colour excess of Tr 16-1 12 is much lower than 
for FO 15, adopting Rv =4.15 for Tr 16-1 12 results in a reduction 
of the distance modulus of 20% only. 



5 HD 93343, A PECULIAR BINARY SYSTEM 



Based on 10 low-dispersion spectra, [Levato et ar] ( [I99l] > reported 
RV variations of HD 93343 (08 Vn) although they could not pro- 
vide firm evidence for the multiplicity of the star. Later on, an ab- 
stract reporting a very preliminary SBl orbital soluti on (Pprb ^ 
44days,e~ 0.39, Ki ~ 71kms-i) was presented by |Solivella & 
|Niemela|(T998) . 

Here, we analyse six medium-resolution Boiler & Chivens and 
nine high-resolution FEROS spectra of HD 93343 spread over six 
years. These data reveal the SB2 nature of this star: one of the 
components of HD 93343 has narrow absorption lines (FWHM ~ 
125 km s^^), whilst the companion displays much broader lines 
(FWHM ~ 500 km s^^ ! ). As can be seen on Fig.js] the sharper lines 
obviously move in radial velocity. The RV changes of the broad 
lines are less obvious to see. However simultaneous fits of two 
Gaussians to the observed profiles clearly indicate that the broad 
lines move in phase opposition with the narrow lines, although 
their RVs are affected by large uncertainties (at least 20kms~^) 
due to their widths. Our RV data allow us to establish a prelimi- 
nary estimate of the mass-ratio of the system rrinarrow/w-broad ~ 
0.63 ±0.06. 
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Table 5. Radial velocities of HD 93343 as inferred from a simultaneous fit 
of two Gaussians to the Hel AA4471, 4921. 5876 and Hell A4542 line 
profiles observed in our data. Note that for the Boiler & Chivens spectra 
(B&C), only the He I A 447 1 and He II A 4542 lines were measured. 



HJD-245000() 


Narrow lines 


Broad lines 


Inst. 


534.743 


-52.5 


101.0 


B&C 


535.749 


-55.8 


87.7 


B&C 


536.716 


-47.9 


63.0 


B&C 


537.731 


-32.2 


53.9 


B&C 


538.715 


-50.9 


47.4 


B&C 


539.717 


-38.3 


42.4 


B&C 


1299.673 


4.8 


-20.6 


F 


1300.677 


0.2 


-27.4 


F 


1301.630 


-9.5 


-19.4 


F 


1302.605 


-22.7 


-6.7 


F 


1304.650 


-30.0 


13.0 


F 


1670.632 


-79.2 


70.5 


F 


1672.617 


-81.3 


60.7 


F 


2783.550 


-80.5 


49.1 


F 


2784.554 


-81.3 


59.3 


F 



The apparent EWs (i.e. with respect to the combined primary 
+ secondary continuum) of the He I A 447 1 line as determined from 
the two Gaussian fits are 0.57 ± 0.02 and 0.23 ± 0.05 A for the 
broad and narrow hne component respectively. For the He II A 4542 
lines, the apparent EWs are 0.41 ± 0.05 and 0.17 ± 0.04 A respec- 
tively. From these numbers, we infer a spectral type 08 for both 
stars (logW^' = 0.14 ± 0.06 for the broad line component and 
log W' — 0.11 ± 0.14 for the star with the sharper lines, see |Conti 
|1973b| and Mathys 1988 1. We note that the spectral type of the com- 
ponent with the sharper lines could actually be in the range 07-8.5. 
The same equivalent widths allow us to estimate that the broad line 
component should be roughly 2.4 times brighter in the optical than 
the narrow line component. Unfortunately, an attempt to disentan- 
gle the spectra of the two stars failed because the RV changes sam- 
pled by our data do not cover the full width of the broad lines. In 
such a situation, the disentangling method cannot recover the spec- 
tra of the components (see |Gonzalez & Levato|2006l l. 

The likely combination of a slow and fast rotator in HD 93343 
is reminiscent of the situation encountered in Plaskett's Star 
(HD47129, see |Linder et al.|2008^ and this system could thus pro- 
vide an interesting test case for binary evolution theories (e.g. [de] 
[Mink et al.|2009 ). Unfortunately, the existing data are not sufficient 
to establish a full orbital solution. From the RV changes observed 
over 6 consecutive nights in March 1997 and May 1999 (see Ta- 
ble[5](, we can guess that the orbital period is likely of order a few 
weeks. A very preliminary Fourier analysis indicates the presence 
of several alias peaks of identical amplitude in the periodogram 
for periods around 40 - 50 days, with the highest one correspond- 
ing to a period of 44.54 days. We note that the latter result is in 
good agreement with the preliminary period reported by |Solivella| 
[&Ni eiTiela (19981. 

We finally analysed an extensive set of 6101 y-band photo- 
metric measurements of HD 93343 taken between June 2004 and 
April 2008 with the same instrumentation as for Tr 16-112. Al- 
though the photometric variations are rather small, the |Heck et al.| 
( |1985j periodogram actually presents its highest peak for a pe- 
riod of 44.15 days with a peak-to-peak amplitude of 0.016 mag. 
Whilst another peak with almost identical amplitude is found in 
the periodogram for a period of 80.3 days, we nevertheless con- 
clude that the orbital period of HD 93343 is indeed likely to be 
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Figure 5. Montage of several of the He I A 447 1 profiles observed in the 
spectrum of HD 93343. The date of the observation is indicated on the right 
in the format HJD-2450000. 
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Figure 6. Photometric variations of HD 93343. The top panel illustrates the 
raw data (A V) folded with a period of 44.15 days. Phase zero was arbitrar- 
ily set to HJD 2450000. The bottom panel yields the same light curve, but 
where the data of each single observing night have now been averaged. The 
error bars on the data yield the dispersion of A y about the mean for each 
observing night. 

around 44 days. We stress that the light curve folded with a period 
of 44. 15 days indicates a single modulation per cycle (see Fig.[6| 
unlike what one would expect for purely ellipsoidal variations as- 
suming the system has a circular orbit. A modelling of these pho- 
tometric variations must however await the results of a more exten- 
sive spectroscopic (and photometric) monitoring campaign of this 
system. 
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6 THE RV CONSTANT STAR HD 93250 



Last, but not least, HD 93250 (03.5, |mibom et al.|2002l l was con- 
sidered to have a constant radial velocity by |Levato et al.|fl991| l. 
This statement was based on 9 low-resolution spectra obtained on 9 
consecutive nights. The 1-ct dispersion of their RV measurements 
was 11.6kms~^. However, some indirect indications of a possi- 
ble binarity were found in the radio and X-ray domains. In fact, 
[Leitherer et al.| ( [T995[ l reported on an Australia Telescope Compact 
Array radio observation of HD 93250. The star was detected only at 
the 8.64 GHz frequency. Assuming the radio flux to be of thermal 
origin and adopting a distance of 2.2 kpc, [Leitherer et al.|jT995^ 
inferred a mass loss rate of log M ~ —4.39 ± 0.15, considerably 
larger than the value de rived from analyses of the optical spectra 
of this sta{^ Therefore, Leitherer et al.| |l995'l suggested that the 
radio emission is essentially non-thermal. Since non-thermal (syn- 
chrotron) radio emission is nowadays commonly associated with a 
wind-wind collision in an early-type binary (see e.g. |De Becker] 
[2007| >, one would thus expect HD 93250 to be a binary system. 
Some hints for binarity come also from the X-ray domain. In- 
deed, Chandra observations revealed that this star has the largest 
logLx/iboi (= —6.5) among the objects studied in this paper 
( [Sanchawala et al.]|2007[ l. Moreover, the star has been frequently 
observed with XMM-Newton in the course of various campaigns on 
the Carina region. In the framework of a global study of the X-ray 
properties of OB stars (Naze 2009, in preparation), we have re- 
trieved and analysed 21 XMM-Newton-EPlC spectra of HD 93250 
from the 2XMM catalogue (Wats on et n^|2009| l. The X -ray flux 
of HD 93250 displays clear variability between the different point- 
ings (see Fig.|7]l, whilst most single 0-type stars usually exhibit a 
rather constant X-ray flux. The EPIC spectra of HD 93250 can be 
fitted with a two-temperature optically thin thermal plasma mekal 
model with temperatures around 0.20 keV and 1 .7 - 2.8 keV. These 
results, i.e. the existence of a relatively hard thermal component and 
the significant X-ray variability could be seen either as further evi- 
dence for HD 93250 being a colliding wind binary or alternatively 
they could indicate that t his star harbours a magnetically confined 
wind similar to 61^ Ori C (|Gagne et al.|2005 i. 



Six FEROS high-resolution spectra of HD 93250 were anal- 
ysed. The spectrum of this star displays several emission lines: 
N III AA 4634-40, C ill A A 5696, 672 1 , 6727-3 1 and Si I v A A 4089, 
4116 are the most prominent ones. A weak Niv A 4058 emission 
is also seen, whilst He II A 4686 is in strong absorption (EW — 
0.54 ± 0.06 A) with a slightly asymmetric profile. These features 
indicate that HD 93250 is an O V((f+)) star. |Walbom et al.|j2002| 
introduced the new 03.5 spectral type and proposed HD 93250 to 
be a prototype of this category. This result is in good agreement 
with the ratio of the equivalent widths of the He I A 4471 (EW — 
0.119 ± O.OOSA) and Hell A 4542 (E W = 0.690 ± O.OISA) 
classification lines. The |Conti i 1973b^ criterion yields log W' — 



0.76 ± 0.03 corresponding to an 04 spectral typ^ In summary, 
HD 93250 hence appears to be an 03.5 V((f+)) star. 

Our RV measurements indicate l-cr RV dispersions of less 
than 2.0kms-i for the Hel A4471, Hell AA4542 and 5412 ab- 
sorption lines. The same conclusion holds if we consider the mean 
RV of seven lines (a = 1.3kms^^, see Table|6]l. Note that the 
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Figure 7. From top to bottom the various panels show the observed X-ray 
flux of HD 93250 (in units 10" ergcm'^ g-i) over the 0.5 - lOkeV do- 
main, the wind hydrogen column density (in units 10^^ cm~^) in addition 
to the interstellar column of 2.55 X 10^^ cm~^, the temperatures of the first 
and second thermal component (in keV). 



Table 6. Radial velocities of HD 93250 as measured from the He I A A 447 1 , 
5876, Hell AA4026, 4200, 4542, 5412 and OlII A5592 absorption lines. 



HJD-2450000 Mean RV 



1304.551 


-9.9 


2782.538 


-9.1 


2784.526 


-12.5 


3130.501 


-9.2 


3131.455 


-9.3 


3133.454 


-10.5 



largest variability (cr = 6.5kms~^) is found for the Hell A 4686 
line. This is no surprise since this line shows some hints for a P- 
Cygni type profile and the variability is thus most likely due to 
variations of the stellar wind. For comparison, the RV dispersion 
of the diffuse interstellar band at 5780 A amounts to 0.7 km s"^. 
Therefore, we conclude that no significant RV variations are ob- 
served for the photospheric lines of HD 93250. Our measurements 
sample time scales of 3 - 5 days as well as of several years. The 
lack of significant RV variations therefore indicates that HD 93250 
is most likely not a spectroscopic binary with an orbital period of 
order a few days to a few weeks or a few years. Finally, it should be 
noted that |Mason et al.|(l99 8) did not report a visual companion to 
HD 93250. In conclusion, HD 93250 is probably either a single star 
or a moderately wide binary system seen under a low inclination 
angle. 



S Indeed, 



Martins et al. 



2005b 



subsequently derived log KI = —6.25 it 



0.7 from a model atmosphere fit. 

^ Note that Conti ( 1973b^ did not consider an 03.5 spectral type in his 
classification scheme. 



7 O-TYPE BINARIES IN AND AROUND THE 
TRUMPLER 16 CLUSTER 

In this section, we briefly discuss some of the results obtained so 
far in the course of the optical campaign of the X-Mega project. 
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In addition to the stars studied in this paper, we include the re- 
sults on some other stars: Tr 16-34, HD 93205, Tr 16-104, Trl6- 
1 10, HD 93 161, FO 15 taken from the previous papers of this series 
I Albacete Colombo et al.|2001|[2002l|Morrell et al.|2001[[Rauwetl 



|al.|200H|Naze et al.|2005[|Niemela et al.|2006) as well as Tr 16 
1 and HD 93403 from the work by |Freyhanuner et al^ j2001) and 
|Rauw et ar] ( |2000l > respectively. 

The properties of the eclipsing systems, Trl6-1, Tr 16-104 
and FO 15 (Freyhammer et al.|[200T| |Rauw et al.][200H |Niemela| 
|et al. [2006) 1, are particularly useful since their light curve and or- 
bital solution provide absolute parameters and allow us to constrain 
their distance. Indeed, the investigations of these three eclipsing 
systems agree with a distance of about 2.5 kpc. For the ellipsoidal 
variables, Tr 16-112 and HD 93205, the uncertainties on the orbital 
inclination and hence on the absolute parameters are much larger 
( |Antokhina et al.|2000| and this work). 

For the non-eclipsing binaries and for the probably single star 
HD 93250, we have assumed a distance of 2.5 kpc. For each star, we 
have evaluated the effective temperature and the bolometric cor- 
rections according to [Martins et al.| j2005a| l and [Martins & Plez| 
( |2006) . We have generally assumed an uncertainty of half a spec- 
tral type. To derive the luminosities, we need the observed mag- 
nitude, the reddening, the brightness ratios between the different 
components, the bolometric corrections and an estimate of the dis- 
tance. The observed magnitudes and colours were taken from the 
WEBDA database. The relative brightness of the components of the 
binary systems was estimated in the various papers from the ratio 
of the EWs of the primary and secondary lines. The reddening was 
evaluated according to the results of [Carraro et al.[ ( [20041 . These 
authors argued that a unique reddening law is not appropriate to 
study the whole Carina region and they estimated different values 
of the selective extinction Ry = Ay /E{B — V) for each cluster. 

This approach yields the bolometric luminosities quoted in Ta- 
ble[8] For Tr 16-1 12, we use the luminosities derived from the light 
curve analysis, rather than assuming a distance of 2.5 kpc. The lu- 
minosity of HD93161B might be overestimated if the star turned 
out to be indeed another binary system (see [Naze et ^L|[2005| >. 
Note that for HD 93403, the evolutionary masses are revised down- 
wards compared to the previous analysis ( [Rauw et al.[[2000] > and 
the evolutionary mass ratio is now in better agreement with the ob- 
served value (1.75 ± 0.06). Finally, we stress that the bolomet- 
ric luminosity of HD 93250 that we obtain through our approach 
(log L/L© = 6.09 ± 0.07) is in very good agreement with the lu 
minosity logL/Le — 6.12'^'^^ inferred by 



Martins et al. 



(2005b 



using a different method. 

In Fig.jS] we have plotted the various stars of TablejS] in a 
Hertzsprung-Russell diagram along with the evolutionary tracks 
from [Meynet & Maeder| ( [2003[ l for solar metallicity. Provided 
that all stars (except maybe Tr 16-112) are indeed at the same 
distance from Earth, this diagram suggests that the components 
of HD 93161 and HD 93403 are evolved off the zero age main- 
sequence with an age of about 2-3 Myr, unlike the components 
of the eclipsing binaries and of Tr 16-110 which are found very 
close to the zero-age main-sequence (ZAMS). 

Using all the available orbital solutions for early-type binary 
systems in or around Trumpler 16, we constructed the period - ec- 
centricity diagram shown in Fig.[9] This figure confirms our asser- 
tion from Fig.[8]that HD 93 161 A must be older and more evolved 
than the main-sequence stars in the core of Tr 1 6. In fact, for the lat- 
ter, significant eccentricities are found for orbital periods as short 
as ~ 4 days, indicating that tidal interactions have not yet had the 
time to circularize the eccentric orbits. 



Table 8. Revised temperatures and bolometric luminosities of the compo- 
nents of eclipsing binaries (top part), the ellipsoidal variables (second part), 
the non-eclipsing binaries (third part) and two presumably single 0-stars 
(lower part) in the Carina complex. 



System 


Spectral type 


Tcfl (K) 


los L /Lo 


Tr 16-la 


09.5V 


32000 ± 500 


4.42 ±0.06 


Tr 16-lb 


B0.3V 


30000 ± 800 


4.14 ± 0.07 


Tr 16-104a 


07V 


36900 ± 1000 


4.92 ±0.11 


Tr 16-104C 


09.5V 


31900 ± 1000 


4.34 ± 0.09 


F015a 


05.5 V 


40000 ± 1000 


5.09± 0.01 


F015b 


09.5 V 


32000 ± 500 


4.41± 0.01 


HD 93205a 


03.5 V 


43900 + 1000 


5 32 + 1 4 


HD 93205b 


08V 


34900 ± 1000 


4.71 ± 0.10 


Trl6-112a 


05.5-6V 


39000 ± 1000 


5.57± 0.04 


Trl6-112b 


B2V-III 


22000 ± 1500 


3.83± 0.16 


HD93161Aa 


08 V 


34900 ± 1000 


5.16 ±0.04 


HD93161Ab 


09V 


32900 ± 1000 


4.85 ±0.06 


HD 93343a 


08V 


34900 ± 1000 


4.86 ±0.12 


HD 93343b 


07-8.5V 


34900 ± 1000 


4.48 ±0.12 


HD 93403a 


05.5I 


37700 ± 900 


5.72 ±0.05 


HD 93403b 


07 V 


36900 ± 1000 


4.97 ± 0.22 


Tr 16-34a 


07.5 V 


35900 ± 1000 


4.89 ±0.09 


Tr 16-34b 


09 V 


32900 ± 1000 


4.44 ± 0.12 


Trl6-110a 


07 V 


36900 ± 1000 


4.86 ±0.06 


Trl6-110b 


08 V 


34900 ± 1000 


4.69 ±0.06 


Trl6-110c 


09V 


32900 ± 1000 


4.55 ±0.07 


HD93161B 


06.5V(f) 


37900 ± 1000 


5.50 ± 0.04 


HD 93250 


03.5V((f+)) 


44000 ± 300 


6.09 ±0.07 



O 



ZAMS 




4.8 



Figure 8. Location in the Hertzsprung-Russell diagram of the components 
of various early-type systems of the Carina region. Eclipsing binaries are 
shown by filled dots, whilst ellipsoidal vaiiables are shown as filled trian- 
gles. Non-eclipsing binary systems are indicated by open squares and the 
asterisks stand for two probably single stars. The single star evolutionary 
tracks from [Meynet & Maeder| j2003) for an initial rotational velocity of 
300 km s ~ ^ at solar metallicity are shown for those parts of the tracks where 
the surface hydrogen abundance X 0.4. 
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Table 7. Orbital properties of the short-period spectroscopic binaries in the Carina complex. For the third component of the triple system Tr 16-104, the two 
possible solutions are listed (see |Rauw et al.|2001^ . In the last column, SBl and SB2 stand respectively for a single-lined and a double-lined spectroscopic 
binary, whilst the e and v letters indicate systems that are known to display photometric eclipses or ellipsoidal variations respectively. 



System 


spectral types 


Porb 


e 




q = m2/mi 


Type 






(days) 




n 






Trl6-1 


O9.5V-hB0.3V 


1.4693 


0.0 




0.83 ±0.08 


SB2e 


Tr 16-34 


08V-H09.5V 


2.3000 


0.0 




0.76 ±0.01 


SB2 


HD 93205 


03V-H08V 


6.0803 


0.371 ± 0.005 


50.6 ±0.9 


0.42 ±0.01 


SB2v 


Tr 16 - 104 


O7V + O9.5V + B0.2IV 


f 2.1529 


0.0 




0.65 ±0.03 


SB2e 


287 ± 20, 238 ± 8 


1 285,1340 


0.25 ±0.05, 0.37 ±0.04 




SBl 


Trie - 110 


f 07V + 08V 


3.6284 


0.06 ±0.03 


355 ± 15 


0.96 ±0.04 


SB2 














1 09V 


5.034 


0.09 ±0.04 


209 ± 30 




SBl 


HD93161A 


08V-H09V 


8.566 


0.0 




0.76 ±0.01 


SB2 


F015 


05.5V-h09.5V 


1.1414 


0.0 




0.52 ±0.04 


SB2e 


Tr 16-1 12 


05.5V-HB2V-III 


4.0157 


0.15 ±0.01 


342.8 ± 5.6 


0.29 ±0.01 


SB2v 


HD 93403 


05.51 + 07V 


15.093 


0.23 ±0.02 


22.5 ±4.4 


0.57 ± 0.02 


SB2 



CIJ 

o 
o 



0.4 



0.3 



0,2 



0.1 



10 100 
Porb (days) 
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Figure 9. Eccentricity of the orbits of a sample of early-type binaries in the 
Carina complex as a function of the orbital period. The probable members 
of Tr 16 are shown by filled symbols, whereas the two binary systems out- 
side the cluster (HD 93403 and HD 93161A) are indicated by open symbols. 
For the third component of the triple system Tr 16-104, the two possible so- 
lutions are shown, connected by a dotted line. 



At least three systems of multiplicity three or higher were 
found among the early-type binaries in the Carina complex. Con- 
cerning the fraction of systems of higher multiplicity, |Tokovinin| 
( |2004^ found that about 43% of the nearby low-mass (0.5 - 1.5 Mq) 
spectroscopic binaries with P < 10 days quoted in the Eighth Cata- 
logue of Spectroscopic Binaries have indeed a tertiary component. 
The third component provides a natural sink for the angular mo- 
mentum that needs to be removed from a binary system in order to 
make it a close one. |Tokovinin| ( |2004[ l quotes an empirical stability 
limit for a triple system: 



Pout (1 - Go 



> 5 Pin 



where Pn, Pout and Cout correspond respectively to the orbital pe- 
riod of the close binary, the orbital period of the third star around 
the centre of mass of the close binary and the eccentricity of the 
third star's orbit. Both possible orbital solutions for the third com- 
ponent in Tr 16-104 largely satisfy this criterion, whilst the orbital 
period of the SBl component in Tr 16-110 assumed to correspond 
to Pout of a triple system obviously cannot meet this condition, 
thus confirming the idea that there should exist a fourth (unseen) 
component in this system. 



8 CONCLUSIONS 

In this paper, we have presented the very first SB2 orbital solution 
of the Tr 16-1 12 binary system. We derived a rather large mass ra- 
tio and confirmed that this binary has an orbit with a significant 
eccentricity despite of its short orbital period. This system displays 
ellipsoidal variations that can be explained by the primary filling 
about 73% of its Roche lobe at periastron. Furthermore, we have 
presented the first clear evidence for an SB2 signature in the spec- 
trum of HD 93343. This system likely harbours one fast and one 
slow rotator and could currently be a post-Roche lobe overflow evo- 
lutionary stage, making it a somewhat less extreme version of Plas- 
kett's Star. Finally, we have shown that the properties of the com- 
ponents of most of the binaries in the Carina region (in or around 
Trumpler 16) suggest that these are very young objects where tidal 
interactions have not yet had the time to circularize the orbits. One 
point that remains to be studied in the forthcoming years is the lack 
of a clear binary signature in HD 93250. This system displays a 
non-thermal radio emission and such a phenomenon is a priori not 
expected for a single star. Interferometric observations could help 
to search for a companion that might be either too faint or too far 
away from the 03.5 primary to be detected in spectroscopy. 
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